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Abstract 


In  the  case  of  the  zero  energy  scattering  of  one  compound  system  by 
another,  where  one  real  scattering  length  completely  characterizes  the  problem, 
(e.g.,  the  reaction  A+B  ->  C+D,  in  addition  to  A+B  ->  A+B,  cannot  take  place) 
it  has  previously  been  shown  that  the  Kohn-Hulthen  variational  principle 
provides  an  upper  bound  on  the  scattering  length  if  no  composite  bound  states 
exist.  The  extension  of  this  result  to  the  case  where  one  or  more  composite 
bound  states  do  exist  is  presented  here.  The  inclusion  of  tensor  forces, 
exchange  forces,  and  Coulomb  forces  is  allowed.  Several  methods  are  given 
for  obtaining  a  rigorous  upper  bound  on  the  scattering  length,  which  involve 
the  addition  of  certain  positive  terms  to  the  Kohn-Hulthen  variational  ex- 
pression. The  approximate  information  about  the  composite  bound  states  which 
is  required  to  construct  these  additional  terms  can  be  found  by  standard  methods. 
As  a  consequence  of  one  of  the  results  obtained,  it  is  shown  that  under  certain 
circumstances  some  ordinary  variational  calculations  give  a  bound.   Thus,  an 
analysis  of  a  previous  calculation  in  the  light  of  the  present  results  leads, 
without  further  calculations,  to  a  rigorous  upper  bound  on  the  singlet  electron- 
hydrogen  scattering  length. 
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I.  INTRODUCTION 
The  value  of  variational  methods  in  scattering  problems  is  greatly  increased 

when  these  methods  are  based  on  a  minimum  (or  maximum)  principle.  This  is  especially 
true  for  the  difficult  case  of  scattering  by  a  compound  system,  where  calculations 
in  which  different  trial  functions  are  used  can  lead  to  quite  different  results. 
Such  a  minimum  principle,  for  zero  energy  scattering,  has  already  been  obtained  for 
the  case  in  which  the  scattered  particle  cannot  be  bound  to  the  scattering  system  . 
The  present  paper  is  concerned  with  the  extension  of  this  result  to  the  case  in 
which  one  or  more  composite  bound  states  exist.  For  clarity  of  presentation  the 
detailed  discussion  is  confined  to  the  problem  of  the  scattering  of  a  spinless, 
neutral  particle,  of  zero  orbital  angular  momentum,  from  a  short  range  center  of 
force.  The  generalization  to  the  case  where  the  scattering  is  by  a  compound  system, 

and  where  long  range  Coulomb  forces  exist,  is  identical  with  the  corresponding 

2 
generalization  in  the  problem  in  which  no  composite  bound  states  exist  ,  and  the 

details  will  be  omitted.  Further,  the  inclusion  of  tensor  forces  is  allowed. 

The  formulation,  given  in  II,  for  the  scattering  of  a  particle  by  a  compound  system 

can  be  extended,  in  a  straightforward  way,  to  treat  the  scattering  of  one  compound 

system  by  another.  It  will  in  fact  be  understood  in  the  following  that  unless 

otherwise  stated  (in  particular,  the  result  of  Section  IV  is  excepted)  each  of 

the  results  obtained  has,  under  the  condition  discussed  in  II,  direct  applicability 

to  this  wider  class  of  problems  (the  essential  requirement  is  that  the  scattering 

be  completely  characterized  by  one  real  scattering  length).  The  problems  which 

may  be  treated  include  zero  energy  scattering  of  a  neutron  or  a  proton  (or  in 

fact  of  nuclei)  by  nuclei,  assuming  realistic  nuclear  potentials.  Similarly,  the 

scattering  of  electrons (or  atoms^by  atoms  may  be  treated. 

As  discussed  in  II,  the  scattering  cannot  be  characterized  by  a  real 
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scattering  length  if  the  radiative  capture  process  is  possible.  We  shall,  in 
the  following,  make  the  approximation  of  ignoring  the  interaction  of  the 
particle  with  the  radiation  field.  The  bound  on  the  scattering  length  is  then 

rigorous  only  to  the  extent  that  this  generally  excellent  approximation  is 

3 
in  fact  valid. 

In  Section  II  a  number  of  results  are  obtained  which  provide  upper  bounds 
on  the  scattering  length  for  the  case  where  only  one  bound  state  exists.  In 
Section  III  it  is  shown  that  there  are  well  defined  circumstances  for  which  some 
ordinary  variational  calculations  provide  a  bound.  An  alternate  method  for 
obtaining  bounds,  based  on  the  Kato  formalism,  is  considered  in  Section  IV.  The 
extension  to  the  case  of  many  bound  states  is  treated  in  Section  V.  Since  in 
this  case  the  exact  scattering  function  has  a  number  of  nodes  (even  for  the 
scattering  of  a  particle  by  a  center  of  force),  making  the  construction  of  an 
accurate  trial  function  more  difficult,  it  would  seem  that  a  minimum  principle 
for  the  scattering  length  is  even  more  valuable  here. 

II.  THE  CASE  OF  ONE  BOUND  STATE 

We  begin  by  considering  the  problem  of  the  zero  energy,  zero  orbital  angular 
momentum  scattering  of  a  spinless  neutral  particle  of  mass  m  by  a  short  range 
center  of  force  whose  strength  is  such  that  one,  and  only  one,  bound  state  exists. 
As  shown  in  I,  the  scattering  length  A  may,  in  general,  be  written 

A  -  A.  -  I  u.Zu.dr  +  fw^wdr  ,  (2.1) 

where  u.  is  a  trial  scattering  function,  w  is  the  difference  between  u.  and 

the  exact  function,  u,  and  <£  is  defined,  in  terms  of  the  Hamiltonian,  H,  as 

2 
X   -  -(2mA2)H  -  2—  -  (2raA2)v(r)  . 
dr 

(Actually,   since  we  write  the  one  dimensional  equations  with  a  factor  of  l/r 
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separated  out  from  the  wave  function,  the  differential  operator  in  H  is  not 
exactly  the  kinetic  energy  operator.)  In  Eq.  (2.1),  and  in  the  following,  all 
one  dimensional  integrals  are  understood  to  have  the  limits  zero  and  infinity. 
The  exact  function  is  determined  by  the  equation 

X   u  =  0  , 
with  boundary  conditions 

u(0)  -  0  , 

/  \  *  (2.2) 

u(r)  ->  A  -  r,    for  r  ->  oo  • 

The  trial  function,  u. ,  satisfies  the  same  boundary  conditions,  with  A  replaced 

W 

by  A.  .     It  was  shown  in  I  that  if  no  bound  state  exists,   and  if  the  normalization 
given  by  Eqs.    (2.2)    (referred  to  in  the  following  as  the  " appropriate" 
normalization)  is  employed,  then 

f  w<£w      <     0  ,  (2.3) 

so  that  the  Kohn-Hulthen  variational  expression  provides  an  upper  bound  on  A. 
The  inequality,  Eq.  (2.3)  follows  from  the  fact  that  under  the  circumstance 
considered,  the  Hamiltonian  is  positive  definite  on  a  space  of  quadratically 
integrable  functions  or,  more  generally,  (see  I)  for  functions  which  approach 
some  finite,  possibly  non-vanishing  constant  at  infinity.  Thus  the  inequality 
will  not  generally  be  valid  if  the  asymptotic  form  f or  u  has  >*iat  we  shall  refer 
to  as  the  "  inappropriate"  normalization,  namely, 

u  ->  1  -  r/A,  for  r  ->  oo.  (2.U) 

We  assume  in  the  following  that  the  "appropriate"  normalization  is  employed,, 
(It  is  of  course  true  that  at  non-aero  energies  one  cannot  generally  say  which 
normalization  is  superior,  and  that  even  at  zero  energy  the  "  inappropriate" 
normalization  provides  a  variational  estimate  for  if  not  a  bound  on  the 
scattering  length.) 
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In  the  case  where  one  bound  state  exists,  which  is  our  present  concern,  the 
inequality,  Eq.(2.3),  will  not  be  valid  unless  the  error  function,  w,  is  orthogonal 
to  the  exact  bound  state  function,  denoted  by  u  .,  where  (u  nu  ,)  =1.  [We  have 
oCu,  ■  e.  u  .,  where  e..  is  related  to  the  bound  state  energy  E  by  e."  -  (2m/h  )E^j 
Since  the  exact  solutions  u  and  u  .  are  orthogonal,  the  difference  function 
w  will  be  orthogonal  to  u  .  provided  u.  and  u  -  are  orthogonal.  We  therefore 
have  the  result  that  the  Kohn-Hulthen  variational  principle  gives  an  upper 
bound  on  the  scattering  length  provided  a  trial  function  is  used  which  is 
orthogonal  to  the  exact  bound  state  wave  function* 

It  is  easily  shown  that  this  statement  is  equivalent  to  the  expression 
for  the  bound  given  below  in  Eq.  (2.5),  where  the  only  requirement  on  u.  is  that 
it  satisfy  the  correct  boundary  conditions.  However,  we  will  actually  obtain 
Eq.  (2.5)  in  a  slightly  different  way,  which  is  the  more  natural  generalization 
of  the  method  used  for  the  no  bound  state  case.  Although  w  is  in  fact  unknowi, 
we  may  formally  construct  a  function,  namely, 

[w  -  (w,uel)u£l], 
which  is  orthogonal  ton  ..  It  then  follows,  exactly  as  in  the  case  for  no 
bound  states,  that 

([w-(w,u6l)u6l]/  fr-(w,»A>»J)  <  °- 

This  inequality  may  be  written 

(w,£w)  <  e1(ut,uel)2^ 

where  we  have  made  use  of  the  relation  (w,n  ..)  ■  (u.,u  ..).  Thus,  an  upper 
bound  on  A  may  be  obtained  from 

A  ±  At  "  /  V  ut  dr  +  «1(J  utueldr)2'  &.*) 
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The  difficulty  with  Eq.  (2.5)  is  that  u   is  not  known  in  general  . 
Therefore  unless  the  last  term  in  Eq.  (2.5)  is  small,  for  a  given  (presumably) 
reasonably  accurate  trial  function  u,  ,  and  shows  very  small  variation  as  the 
accuracy  of  the  trial  bound  state  function,  u  , . ,  is  increased,  one  cannot 
claim  to  have  a  bound.  In  fact,  preliminary  calculations  based  on  Eq.  (2.5), 
for  the  singlet  e~H  problem,  gave  negative  results  in  the  sense  that  the  last 
term  in  Eq.  (2.5)  was  large  and  showed  large  variations  for  a  series  of  functions 
u  ...  ranging  from  a  two  parameter  to  an  eleven  parameter  function  with  u.  that 
given  by  Seaton.   The  failure  of  this  method  in  the  particularly  favorable 
e~H  case  indicates  that  Eq.  (2.5)  is  not  a  useful  form.  However,  from 
Eq.  (2.5)  one  can  obtain  an  expression  which  does  not  depend  on  the  unknown 
function  u  .  by  writing 

htf  Vadr'2  "  \  </  Vueidr)2 

where  use  has  been  made  of  the  Schwarz  inequality.     We  then  obtain 

A  <  At  -  JutAtdr  +  i-  J(£ut)2dr  .  (2.6) 

The  experimental  value  of  the  binding  energy,  if  available,  gives  e,. 
Alternatively,  e.  may  be  replaced  by 

"it  "/Vt^elt^V 

where  u  ,.  is  a  normalized  trial  bound  state  function,  while  still  preserving 
elx 

the  inequality,  Eq.  (2.6).  For  the  case  of  scattering  by  a  center  of 
force,  where  quite  accurate  trial  scattering  functions  may  be  constructed 
without  great  difficulty,  Eq.  (2.6)  may  be  expected  to  give  useful  upper 
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bounds  on  A.  rlowever,  the  use  of  the  Schwarz  inequality  in  passing 

from  Eq.(2.50  to  Eq.  (2.6)  leads  to  a  very  great  overestimate  of  the  error 

term  for  any  but  an  excellent  u.  ,     Thus,  the  trial  function  given  by 

Seaton  for  the  singlet  e  H  problem  ,  which  is  expected  to  be  reasonably  accurate, 

leads  nevertheless,  when  used  in  the  appropriate  generalization  of  Eq.  (2.6) 

to  scattering  by  a  compound  svstem,  to  an  upper  bound  on  the  scattering  length 

which  is  so  large  as  to  be  quite  useless.  It  does  not  seem  that  there 

exist  other  many  body  scattering  problems  for  which  a  similar  treatment  would 

yield  more  favorable  results. 

It  was  therefore  necessary  to  consider  a  somewhat  different  approach  to 
the  problem,  one  leading  to  a  bound  on  A  which  does  not  require  a  knowledge 
of  the  exact  composite  bound  state  function  and  which  does  not  utilize  the 
often  very  crude  Schwarz  inequality.  We  now  show  how  such  a  result  can  be 
obtained  based  on  the  following  two  theorems  which,  for  the  purposes  of 
later  generalization,  are  quoted  in  a  more  general  form  than  is  needed  in 
this  section. 

Theorem  1:  We  consider  a  system  with  exactly  N  bound  states.  If  a  KxK 
matrix  of  the  Hamiltonian  operator  H  is  formed  (K  may  be  larger  than,  equal 
to  or  smaller  than  N)  using  a  set  of  K  orthonormal  functions,  and  is 
diagonalized,  with  the  diagonal  elements  arranged  in  increasing  order,  i.e., 
-(K)   F(K)      F(K)       E(K) 

then  E.   >  E.,  if  i  <  N.  Here  E.  is  the  energy  of  the  i—  bound  state.  For 
i  >  N,  we  have  e|K'  >  0. 

Theorem  2:  If  the  Hamiltonian  matrix  is  constructed  in  the  same  way  as 
discussed  in  Theorem  1,  but  of  rank  K  +  1  with  the  first  K  functions  un- 
altered, then  the  eigenvalues,  E^    ,  satisfy  the  relations 


-  8  - 


E,  <   E<K+1>  <  EfK)  . 
i  —  1    —  1 

We  now  apply  these  theorems  to  the  case  where  only  one  bound  state  exists. 

We  introduce  the  orthonormal  functions  v.  and  v_  and  form  the  2x2  Hamiltonian 

matrix 


where 

H±i  b  (v.,^),       i,J  -1,2. 

It  is  assumed  that  v.  is  a  sufficiently  accurate  ground  state  trial  function 

so  that  /_  -j 

Hll  =  El  K   °# 
(If  this  last  inequality  is  not  satisfied,  the  method  described  below  will 

not  be  applicable.  This, however,  does  not  represent  a  serious  limitation 

since  a  trial  function  with  the  required  accuracy  can  almost  always  be  found.) 

(2)      (2) 
It  then  follows  that  the  eigenvalues  of  the  above  matrix,  E^   and  E^  , 


satisfy  the  relations 


from  Theorem  1,  and 


42)>o, 


E<2)  <  E*1^  0, 


from  Theorem  2,  so  that 

e{2>e2<2)<o. 

Since  the  product  of  the  eigenvalues  of  a  matrix  is  equal  to  the  determinant 
of  the  matrix  we  obtain  the  inequality 

H11H22  ""Li0-  (2-7) 

We  note  that  since  v?  appears  exactly  twice  in  each  of  the  two  terms  of  this 
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inequality,  neither  the  magnitude  nor  the  sign  of  the  normalization  of  v? 
is  relevant,  though  of  course  v?  mu3t  be  normalizable.  (Clearly  the  same 
is  time  of  v  but  it  will  be  convenient  in  the  following  to  think  of  v.  as 

normalized  to  unity.)  We  now  observe  that  the  inequality,  Eq.  (2.7),  is  valid 

i 

if  v?  is  replaced  by  a  normalizable  function  v_  which  is  not  orthogonal  to 

v_ .  This  is  verified  by  writing 

V2  *  V2  "  (VV1)V1  » 
[so  that  (v?,v-)  ■  0  is  satisfied]  and  substituting  this  expression  for  v? 

in  the  integrals  which  appear  in  Eq.  (2.7 )•  Since  the  normalization  of  v_ 

as  well  as  the  value  of  (v?,v..)  are  both  arbitrary,  as  is  clear  from  the  above 

discussion,  we  are  now  in  the  position  to  choose 

v    "  w(X)  s  we    , 
with  X  >  0.  We  then  have  the  inequality  which  is  the  new  form  of  Eq.  (2.7) > 

Jw(X)/w(X)dr  <  JL  (f  u^  w(X)dr)2,  (2.8) 

2   (l) 
with  e_,  ■  -(2m/fo  )E.V   and  u  ...  ■  v..  Since  both  sides  of  the  inequality, 

Eq.  (2.8),  are  continuous  at  X  ■  0  (see  I)  we  have 

/w^wdr  <  (I  u  -..oowdr)     . 
-  elt   V      elt 

Therefore,   since  dfw  ^aLu.,  the  expression  which  provides  an  upper  bound  on  the 
scattering  length  becomes 

A  ±  At  "  /  \l  utdr  +  I^(/  Uelt^  S*>2«  (2-9> 

We  have  then  succeeded  in  our  purpose  in  that  we  have  obtained  a  bound  which 
does  not  require  a  knowledge  of  the  exact  bound  state  function  nor  of   its 
energy,   and  which  has  at  the  same  time  avoided  the  use  of  the  rigorous  but 
often  very  crude  Schwarz  inequality.      (In  the  case  where  the  scatterer  is  a 
compound  system,   the  ground  state  wave  function  of  the  scatterer  must  still 
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n 

be  known  exactly  in  order  to  obtain  a  rigorous  bound.  )  Note  that  the 
earlier  results,  Eqs.  (2. 5)  and  (2.6),  may  be  considered  as  special 
cases  of  Eq.  (2.9);  the  choice  u  ,.  »  u  .  in  Eq.  (2.9)  yields  Eq.  (2.5),  and 
the  use  of  the  Schwarz  inequality  in  the  last  integral  of  Eq.  (2.9)  leads 
back  to  Eq.  (2.6).   It  is  of  course  clear  that  for  a  given  form  of  the  trial 
functions  u.  and  u  ,.  the  optimum  choice  of  the  variational  parameters  is 
such  as  to  minimize  the  right  hand  side  of  Eq.  (2.9),  subject  to  the 
requirement  that  j  u  -,+X  u  -.  +  dr  >  0.  . 

For  the  case  of  arbitrary  values  of  the  orbital  angular  momentum,  L,  the 
inequality  which  corresponds  to  that  given  in  Eq.  (2.9)  has  the  same  form  as 
the  result  for  L  ■  0.  We  merely  replace  A  (with  a  similar  replacement  of  A . ) 
by  Aj.,  defined  in  terms  of  the  phase  shift  17  .  as 

A^   -  -  [lx3x...x(2L+l)]2(tan^L/k2L+1)kii0  . 

Further,  the  trial  function,  u, .  (r),  satisfies  the  boundary  conditions 

uLt(0)  -  0  , 
uu(r)  ->  ALtr"L/(2L+l)  -  rL+1,     for   r  ->  oo, 

n  0 

and^Cnow  contains  the  additional  term  L(L+l)/r  .  Of  course  for  a  particular 
value  of  L  the  bound  state  considered  is  one  with  orbital  angular  momentum  L. 
These  results  for  arbitrary  values  of  L  can  be  extended  in  a  straightforward 
manner  to  the  case  of  scattering  of  one  compound  system  by  another,  where  each 
system  may  carry  a  net  charge.  The  regular  and  irregular  solutions  of  the 
"free"  wave  equation  in  the  limit  of  vanishing  energy  have  previously  been 
given  for  this  general  case  *  The  asymptotic  form  of  the  wave  function  is  taken 
as  the  "appropriate"  linear  combination  of  these  solutions.  The  analysis 

then  proceeds  in  a  manner  similar  to  that  already  described  in  I  and  II  and  we 

o 
shall  omit  further  details  here. 
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III.   CONNECTION  WITH  FORMS  OF  THE  VARIATIONAL  PRINCIPLE 
The  inequality,  Eq.    (2.9),   has  the  particular  consequence  that  there  exist 
well  defined  circumstances,   which  will  now  be  described,   under  which  an  upper 
bound  on  the  scattering  length  may  be  obtained  from  an  ordinary  variational 
calculation.      (As  in  Section  II  the  discussion  in  this  section  will  be  confined 
to  the  case  where  only  one  bound  state  exists.)     We  consider  first  the   Kohn-Hulthen 
principle,  which  has  the  zero  energy  form 

A~    At  -  JutAtdr,  (3.1) 

where  u.    satisfies  the  boundary  conditions  given  in  Eqs.    (2.2).     Now  suppose 
u.    may  be  written, 

ut  =  ut  +buelt  ,  .  <3'?) 

where  u  ...    is  a  trial  bound  state  function  which  is  sufficiently  accurate  to 
give  binding,   i.e., 

elt  =/uelt^uelt  dr>0- 

i 
Since  u  _.  vanishes  at  the  origin  and  at  infinity  u.  satisfies  the  same  boundary 

condition  as  u,  and  can  therefore  itself  serve  as  a  proper  trial  scattering 

function.  Substituting  this  form  for  u.  into  the  variational  expression,  Eq.  (3.l), 

i 

and  determining  b  variationally,   with  u.    and  u  ..    momentarily  considered  as 

fixed,   we  find 

b  ■  -  - —  J  u  ■,  +  *  u.    dr  , 
elt/      elt       t 


and 


A*  At  "  /  ut/ut  dr  +  \rr  (/  uei/  V  dr)2-  (3-3) 


it 

It  follows  from  Eq.  (2.9)  that  the  variational  estimate  of  A  given  bv  Bq.  (3.3) 
is  actually  an  upper  bound.  Conversely,  the  approximation  to  the  zero  energy 
scattering  function,  obtained  in  the  course  of  a  calculation  of  an  upper  bound 
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on  A  based  on  Eq.  (2.9),  should  be  taken  to  be 


\  '   (  Z7  /  uelt^ut  dr)u< 


"It  '   •*«   "     elt 
rather  than  V,  since  it  is  the  former  function  which  can  be  interpreted  as  a 

variationally  determined  trial  function.  We  note  that  the  prescription  given 

by  Kohn  to  evaluate  the  variational  parameters  should  be  used  rather  than  that 

of  Hulthen  since,  while  both  methods  will  provide  an  upper  bound,  for  a 

given  form  of  the  trial  function  the  Kohn  method  will  yield  a  lower  (and 

therefore  better)  estimate  of  the  scattering  length.  (  See  I.) 

It  is  of  interest  to  see  if  there  are  any  variational  calculations  of 

scattering  lengths  reported  in  the  literature  which  were  performed  with  trial 

functions  which  may  be  written  in  the  form  given  by  Eq.  (3.2-'.  Such  calculations 

may  be  reinterpreted  in  the  light  of  the  present  discussion.  As  an  example 

we  note  that  Borowitz  and  Greenberg   have  recently  performed  a  variational 

calculation  for  electron-hydrogen  scattering,  using  as  a  zero  energy  trial 

function  one  which  was  explicitly  constructed  to  be  of  the  form  of  a  trial 

scattering  function  plus  a  multiple  of  a  trial  H~  bound  state  function.  Since 

the  normalization  used  by  these  authors  corresponds  to  what  we  have  called 

the  inappropriate  normalization  [see  Eq.  (2.U)]  neither  the  results  for  the 

singlet  case,  where  one  bound  state  of  H   exists,  nor  for  the  triplet  case, 

where  no  bound  states  exist,  is  necessarily  a  bound.  However,  as  described 

in  I,  for  those  calculations  which  have  been  performed  using  the  inappropriate 

normalization  the  results  may  be  converted,  with  a  trivial  amount  of  labor, 

so  that  they c orrespond  to  the  Kohn-Hulthen  form  and  therefore  do  give  a 

bound  under  the  circumstances  considered.  The  method  of  conversion  may  be 

briefly  restated  as  follows.  If  A.  and  I  =  J  u.cif u.  dr  have  been  evaluated, 
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with  u.    normalized  as  in  Eq.    (2.U),  then  the  bound  on  A  is  given  by 

2 

A  <  A.    -  A     I, 

If  A.  is  positive  and  j A  1 1  is  small  compared  to  unity  then  this  conversion  may 
be  expected  to  give,  in  addition  to  a  bound,  an  improved  approximation  to  the 
scattering  length. 

Returning  to  the  work  of  Borowitz  and  Greenberg,  we  note  that  two  separate 

calculations  of  the  singlet  scattering  length   were  performed,  using  two 

12 

trial  H  functions  obtained  by  Chandrasekhar  •  One  of  these  was  a  two  parameter 

function  containing  no  dependence  on  the  interelectronic  distance,  r  ?,  so  that 
no  polarization  was  allowed  for.  The  other  function,  which  included  a  linear 
r12  dependence,  contained  three  parameters.  In  each  case  the  parameters  used 

where  those  given  by  Chandrasekhar.  The  two  and  three  parameter  functions 

—2  —2  12 

give  values  of  e1+  equal  to  0.0266  a    and  0.0518  a    respectively 
it  o  o 

(a  is  the  Bohr  radius).  The  rigorous  uppper  bounds  obtained  by  the  conversion 
process  described  above  yielded, in  both  cases,  slightly  improved  approximations 
to  the  scattering  length  (see  Table  I). 

IV.  AN  ALTERNATE  METHOD 

In  this  section  we  consider  an  alternate  method  for  obtaining  upper  bounds 
on  the  scattering  length  for  the  case  of  the  scattering  of  a  particle  by  a 
center  of  force  where  one  bound  state  exists.  The  method  is  based  on  the 
formalism  given  by  Kato   for  finding  upper  and  lower  bounds  on  the  scattering 
phase  shift.   (Since  the  Kato  method,  as  opposed  to  the  Kato  identity,  is 
unrelated  to  the  other  methods  discussed  in  this  paper,  this  section  may  be 
read  independently  of  the  others.) 

The  Kato  identity  for  arbitrary  scattering  energy, (hkj/2m,  of  which 
Eq.  (2.1)  is  the  zero  energy  form  with  the  normalization  0  <  9  <  n  (i.e., 
9  /  0),  is 
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k  cot(^-  0)  -  k  cot   (^t-  9)  -  Ju^c/u^dr  +  j  wg^wQdr,      (U.l) 

where  now 

<£.     d      .   (2mA2)  V(r)  ♦  k2  . 
dr 

The  exact  solution,  ufl,   satisfies  the  differential  equation  Xup  *  °>  ancl 


the  boundary  conditions 


tta(0)  -  0  , 


A  -  (U*2) 

ue(r)  ->  sin(kr+7  )/sin(7-  0)  ,     f or  r  ->  oo   , 
where  A  is  the  exact  phase  shift  and  where  9  satisfies  0  <  9  <  n.     The  trial 
function,  u^.  ,   satisfies  boundary  conditions  of  the  form  given  in  Eqs.    (6.2), 
but  with  A  replaced  by  a  trial  phase  shift,  v7  .  ,  which  is  arbitrary;  w_  is 
defined  as  Mq.    -  ufl.     The  Kato  inequalities  are 

"  °  0"1  J  (^uet)2  p"1(ir  -j  V^V1,  ±  P9"1 /(^u9t)2p"ldr#    (U,3) 

Here  p(r)  is  a  non-negative  function  which  must  fall  off  faster  than  l/r 
but  which  is  otherwise  arbitrary  and  is  chosen  for  convenience.     The  numbers 
aQ  and  pfl  are  defined  by  the  associated  eigenvalue  problem,  in  which  the 
equation 

(/+Hp)0Q(r)  -  0 
is  considered.     The  asymptotic  form  of  0e(r)  is  characterized  by  a  phase  shift, 
6(n.).     The  infinite,   discreet  set  of  eignevalues,    j  Ugl   ,  is  defined  such  that 
for  any  integer,   n, 

6(0.^)  -  9  +  nn  . 
The  smallest  positive  eigenvalue  is  aft  and  the  largest  negative  eigenvalue  is  -PQ. 

For  zero  energy  scattering  the  phase  shift  in  the  associated  eignevalue 
problem,   6(m.),  is  given,   according  to  a  theorem  due  to  Levinson,       by 
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6 ((J.)  ■  mn,    if  there  is  no  state  of  exactly  zero  energy,  where  ra  is  the 
number  of  bound  states  for  the  operator  *  +  up.      [if  there  are  ra  bound  states 
and  there  is  one  at  zero  energy  in  addition,  the  phase  shift  is   (m+  xOn.J     we 
are  considering  the  case  for  which  there  is  one  and  only  one  bound  state  for 
the  actual  physical  system,  i.e.,   for  m-  *  0.     We  choose  9  ■  y>  where 

0  <  y  <  ii. (The  point  is  that  we  wish  to  exclude  9-0;   note  that  the  eigen- 
values |i __.  do  not  form  a  discreet  set  at  zero  energy.)     It  is  then  clear  that 

nu 

-p     is  determined  by  the  condition  that  the  operator  at  -p   o  has  associated  with 
it  a  single,   zero  energy  bound  state.     Actually,   according  to  Eq.    (U.3), 
only  a  lower  bound  on  pfl  is  required  in  order  to  obtain  an  upper  bound  on 

1  Wg/w-dr.     Now  a  lower  bound  on  p     is  easily  determined,   according  to  the 
Rayleigh-Ritz  principle,  by 

where  f .    is  an  arbitrary  normalized  trial  bound  state  function.     We  have 

This  result  is  useful  only  if  p       >  0,  that  is,  V.    must  be  sufficiently 

T  *> 

accurate  such  that   (Y.,/^)  >  0. 

If  we  take  the  zero  energy  form  of  Eq.    (U.l),  with  9  ■  y>   and  make  use  of 
Eqs.    (U.3)  and   (U.U),  we  obtain 

A  <  At  -  j  ut/utdr  +   (P^)"1  j  (/ut)2  p_1dr,  (U.5) 

where 

A  «  -   (tany/k)^, 
At  =  -  (tan^/k)^, 

Ut  B   (\*/k)k-0 
and  where  p       is  defined  in  Eq.    (U.h). 
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A  number  of  methods  have  been  given  previously  Ji   J   for  determining  lower 

bounds  on  afl  and  jBQ,  for  arbitrary  scattering  energy.  The  method  just 

described,  while  limited  to  zero  energy  scattering,  has  the  virtue  that  the 

prescription  for  finding  a  bound  on  f3Q  is  quite  general.  The  potential  need 

satisfy  no  special  requirements  other  than  that  it  lead  to  one  bound  state  , 

and  that  a  phase  shift  be  defined.  This  generality  becomes  quite  significant 

in  the  case  of  scattering  by  a  compound  system,  where  methods  based  on  special 

properties  of  the  potential  (such  as  being  everywhere  negative,  for  example) 

will  usually  be  inapplicable.  However,  serious  problems  arise  in  the 

extension  of  Levinson's  theorem  to  the  case  of  scattering  by  a  compound  system, 

and  as  discussed  briefly  in  I,  the  inclusion  of  the  Pauli  principle  raises 

additional  difficulties.  Therefore,  until  these  questions  are  clarified  the 

method  presented  in  this  section  should  be  used  only  for  one  body  scattering 

problems . 

It  is  noted  that  a  ,  which  is  required  to  obtain  a  lower  bound  on  A,   is 
Y 

determined  as  that  value  of  u-  for  which  the  second  bound  state  appears.     Now 

for  a  static  potential,   U(r),the  necessary  condition  for  the  existence  of  n_ 

17 
bound  states,   of  orbital  angular  momentum  L,  is 

(2m/h2)    f  r|0(r)|dr>  ^L+Di^  m  (U.6) 

In  the  case  being  considered,   where  the  actual  potential,  V(r),   supports  one 
bound  state  with  L  =  0,  we  have  from  Eq.    (U.6)  that  a     ,   determined  by 


J  r|-<2in/h2)V(r)  +  a^'pCiOldr  -  2     , 


i 
is  a  lower  bound  on  a  .  This  result  will  be  useful  only  if  a   >  0,  In 

particular,  a  necessary  condition  for  the  applicability  of  the  method  is 

r  V(r)dr  <  2  . 


-(2mA2)  j 
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An  obvious  disadvantage  of  this  approach  is  that    (contrary  to  the  method  for 
determining  p      )   it   provides  no  way  of  successively  improving  the  bound  on 

a  •     Further,   even  if  Levinson's  theorem  can  be  shown  to  hold  in  the  case  of 

Y         ' 

scattering  by  a  compound  system,  no  generalization  of  Eq.  (U.6)  has  been 
obtained;  there  would  still  remain  the  problem  of  determining  conditions  for 
the  non-existence  of  many  body  bound  states. 

V.  THE  CASE  OF  NO RE  THAN  ONE  BOUND  STATE 
The  generalization  of  the  above  results  to  the  case  for  which  there  is 
more  than  one  bound  state  of  zero  angular  momentum  of  the  particle  and  the 
center  of  force  is  fairly  straightforward.  We  assume  that  it  is  known  that 
there  exist  exactly  N  such  bound  states.  The  exact  eigenf unctions  and  the 
exact  energies  of  these  states  will  be  denoted  by  u  1  through  u  „  and  by  E 
through  EN,  respectively,  where 

E1  <  E2  <  ...  EN  <  0  , 
and  e.  =.  (2mAi2)E. ,  1  <  i  <  N.  We  do  not  assume  that  the  eigenf unctions  or 
the  eigenvalues  are  actually  known. 

The  generalization  of  Eq.  (2.6)  is  particularly  simple.  Since  the 

function 

N 
w  sw-F  (w,  u  )u 

1=1 
is  orthogonal  to  each  of  the  bound  states,  it  follows  that 

•  /  i 

w  o  dr  <  0, 

where  it  is  assumed  that  the  "appropriate"  normalization,    namely  that  given  by 

i 

Eq.  (2.2),  has  been  used.  Substituting  for  w  ,  we  find 

kdr)2  . 


/ 


JwXwdr^  i-  (JV^V 
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The  elimination  of  the  unknown  u  .  by  the  Schwarz  inequality  leads  to 

N 


f  wiwdr  <   (  r    -i-)   f  (£uj2dr 
'  i^l  ei      J  % 


However,  an  improved  bound  may  be  obtained  by  writing 

v2 


>', 


i=l 

where  the  summation  in  the  last  term  (which  actually  includes  an  integration 
over  continuum  functions)  is  over  the  complete  set  of  energy  eigenf unctions. 
Making  use  of  the  closure  property  of  this  set  we  have 

j  w/wdr  <j-j  U \)2dr  9  (*.l) 

so  that 

A  <  At  -  j  lyTx^dr  +  ~  f  (/ut)2dr  .  (5.1)' 

If   e     is  unknown,   a  lower  bound  may  be  employed,   obtained,  for  example,  by  the 
method  described  in  Theorem  1. 

An  alternate  proof       of  Eq.    (5.1)  may  be  based  on  the  inequality 

j  1>t(H-c)2  Tt  dr  >  (Ec  -  c)2  j  Vt2  dr  ,  (5.2) 

which  is  verified  by  substituting  the  energy  eigenfunction  expansion  for  the 

normalieable  trial  function,  j*. i  c  *s  an  arbitrary  non-positive  energy  value 

and  E     is  the  energy  eigenvalue  closest  to  c.     The  inequality,  Eq.(5.l), 
c 

follows  from  the  choice  c  ■  E„  /2. 

To  obtain  the  generalization  of  the  more  useful  result,  Eq.    (2.9),   we 
introduce  N  orthonormal  functions  and  construct  and  diagonalize  the  N  x  N 
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matrix  of  the  operator  H,  »as  discussed  in  Theorem  1.  The  N  eigenvalues  of 

CN) 
the  matrix  will  again  be  denoted  by  E  x  ,  while  the  N  orthonormal  eigen- 
vectors of  the  matrix  determine  N  orthonormal  eigenf unctions;  the  latter  will 
be  denoted  by  v. (r).  It  is  assumed  that  the  trial  functions  v. (r)  are 


rW  <  0 


sufficiently  accurate  so  that 

,        1  <  i  <  N.  (5.3) 

We  now  introduce  an  additional  trial  function,  v„  ,(r),  which  is  normalized 

N+± 

and  which  is  orthogonal  to  each  of  the  v. (r),   1  <  i  <  N.     With  this  set  of 
N+l  orthonormal  functions,   we  construct  an   (N+l)  x   (N+l)  matrix  from  the 
operator  H.  Since  the  determinant  of  a  matrix  is  equal  to  the  product  of  its 
eigenvalues,  we  have 


(N) 


«N 


+1,1 


E2(N) 


"n+1,2 


h, 


N+l 


H, 


2,N+1 


.     "n+i, 


N+l 


*&*>*  Bj^x. 


*™.  CM) 


In  the  above  determinant  all  the  elements  not  on  the  diagonal  nor  in  the  last 
row  nor  last  column  are  equal  to  zero,  and  we  havp 

HN+l,i  '  Hi,N+l=  j   \  H  Vl  <*  •  l  ±  1  ±  N+1' 

the  E . ^         ,  1  <  i  <  N+l,   are  the  N+l  ordered  eigenvalues  of  the  matrix. 
Since  we  have  introduced  an  additional  trial  function,   it  follows  from  Theorem 
2  that    E.(N +1'  <  E   ^N  ,   1  <  i  <  N,   and  therefore,   from  Eq.    (5.3),   that 
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E  (N+l)  <   0  fQr  i  <  i  <  n.  On  the  other  hand,  since  there  are  only  N  bound 
states,  it  follows  from  Theorem  1  that  %^+-,  >   °« 


Now  the  expansion  of  the  determinant  in  Eq.  ($.h)   leads  to 

„  (N+l)     v   (N+l) 
^2.  *1 v     N     .  .   (N+l) 

'l  EN 

Since  from  the  previous  discussion  E.^       '/E.         >  0  for  1  <  i  <  N     and 

i     '  l  —   — 


„  (N+l)        (N+l) 
1  \2   1  N  I 

"n+i^+i  =  ^  J^W  (Hi,N+i)  +  E  (NJ   x'"  E  (N)  X  EN+1 


„   (N+lj  >   0^  we  haye  the  inequaiity 


It  is  clear  from  this  last  inequality  that  while  v.,  ,  must  be  normalizable,  its 
specific  normalization  is  irrelevant.  We  choose  for  v„+,  a  function  which, 
as  required,  is  orthogonal  to  the  v.,  1  <  i  <  N,  i.e«, 


v 


N+l 


N 
w  -  ^  (w^v.  . 


(5.6) 


The  difference  function  w  must  be  thought  of  as  normalizable,  but  need  not  be 
orthogonal  to  the  v..  Eq.  (5.5)  then  becomes 

j  w  H  w  dr  >^  -p^  (/  v±  H  ut)2dr,  (5.7) 


Ji 
using  H  w  ■  H  u.,  or,   equivalently, 


(w/wdr<£2    — Tn7  (f  vi^  V^*  (5,7)' 

We  now  obtain  from  this  the  generalization  that  we  want,   namely, 

A  <  At  -  /  u^utdr  +  ^  -^  (J  v/^dr)2   .  (5.8) 

We  might  note  that  Eq.  (5.7)  is  of  precisely  the  same  form  as  Eq.  ($•$), 
with  v„  1  replaced  by  w.  This  invarience  can  perhaps  be  better  understood 
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by  recognizing  that  the  earlier  determinantal  equation,  Eq.  (5.U),  is  invariant 
under  the  same  transformation.  Thus,  the  transformation  from  v.. ,  • . . ,  v„,  v„  - 
to  v_,,..,v„,w,  where  v„  .  is  defined  by  Eq.  (5*6),  is  given  by  the  trans- 
formation matrix  J3,  where 

1 

1 

» 

(wjV^  (w,v2)  '.'l 

The  invariance  of  the  determinantal  equation,  Eq.    (5»U),  under  the  substitution 

vN  i   ->  w  is  then  guaranteed  by  the  fact  that  B^  is  unimodular. 

Finally  we  note  that  the  result  of  Section  III  is  easily  extended  to  the 

case  of  many  bound  states.     That  is,  the  right  hand  side  of  Eq.    (5.8)  may  be 

identified  with  the  Kohn-Hulthen  expression  for  the  s  cattering  length,  in  which 

the  trial  function 

N 
u.    +  Y~    b.v. 
1      fel     i  x 

is  used,  provided  each  of  the  b.  is  determined  variationally. 

VI.  DISCUSSION 
The  two  principle  results  that  have  been  derived,  which  enable  one  to  calculate 
upper  bounds  on  the  scattering  length  for  those  problems  in  which  a  finite 
number  of  composite  bound  states  exist,  are  expressed  by  the  inequalities, 
Eqs.  (5.1)  a«i  (5»8).  The  inequality  of  Eq.  (5.8)  will  lead  to  more  accurate 
estimates  of  the  scattering  length  than  may  be  obtained  from  Eq.  (5*1)  • 
The  latter  result  does  have  one  advantage  in  that  one  has  removed  the  necessity 
of  constructing  a  set  of  trial  bound  state  functions  and  performing  the  various 
integrals  involving  these  functions.  However,  it  may  be  quite  difficult,  in 
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the  more  complicated  problems,  to  construct  a  trial  scattering  function 
sufficiently  accurate  to  obtain  useful  results  from  Eq.    (5.1)  }  the  performance 
of  the  required  integrals  may  also  present  formidable  difficulties.     Eq.    (5.8) 
has  the  interesting  consequence  that  certain  ordinary  variational  calculations 
provide  bounds. 

As  was  mentioned  in  Section  I,  the  results  are  derived  in  terms  of  the 
simple  esse  of  the  zero  orbital  angular  momentum  scattering  of  a  spinless,   neutral 
particle  by  a  center  of  force,  but  may  be  directly  generalized  to  a  far  wider 
class  of  problems.     The  essential  restriction  which  must  be  maintained  is  that 
only  one  exit  channel   (elastic  scattering,  including  Pauli  exchange)  should 
be  open.     As  an  illustration  we  cite  the  extension  of  Eq.    (5.8)  for  the  problem 
of  the  scattering  of  a  particle  by  a  compound  system.     The  appropriate  inequality 
becomes 


*<*t 


-/  It*A£t  iV  *  ^  \  t/t*A  !t  «">2.  («) 


We  have  denoted  the  trial  scattering  function  and  the  trial  function  for  the 

i       composite  bound  state  by  £.   and  J.,   respectively,  instead  of  u.    and  v. 

merely  to  indicate  that  they  represent  the  functions  with  no  factors  of  l/r 

separated  out.     The  operator  TV  is  defined  asA   «  -(2\i/h  )(H-E  ),  where 

S 

H  is  the  total  Hamiltonian,  where  E     is  the  ground  state  binding  energy  of  the 
scattering  system,   and  where  n  is  the  reduced  mass  for  the  scattered  particle 
and  the  compound  scatterer.   If,   asymptotically,  there  exists  a  net  repulsive 
Coulomb  interaction  between  the  scattered  particle  and  the  scattering  system 
then  the  inequality  which  is  applicable  is  identical  in  form  with  that  given 
in  Eq.    (6.1)  but  A  is  now  defined  in  terms  of  the  "additional"    phase  shift,  ^j, 
as20 


A  S-  (tan  *f  /kC2) 


'k-0 
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2 

where  C  is  the  Coulomb  penetration  factor.   (The  asymptotic  form  of  the  trial 

function,  J.  ,  is  altered  in  this  case  but  we  shall  omit  further  details  here.) 
There  is  a  lack  of  completeness  in  the  results  obtained  here  in  the  sense 
that  no  method  has  been  given,  valid  for  many  body  scattering  problems,  which 
provides  a  lower  bound  on  the  scattering  length.  Such  a  method  is  clearly  very 
much  to  be  desired  and  work  is  now  in  progress  on  this  aspect  of  the  problem. 
On  the  other  hand,  while  a  single  bound  on  the  scattering  length  allows  one, 
in  a  particular  calculation,  to  make  no  rigorous  error  estimate  whatsoever,  and 
may  provide  neither  bound  on  the  zero  energy  cross  section,  x.he  knowledge  that 
one  has  'a  minimum  principle  allows  for  successive  improvements  not  only  in  the 
calculated  scattering  length  but  the  trial  wave  function  as  well.  Since  an 
accurate  estimate  of  the  effective  range,  r  ,  (but  not  the  shape  dependent 
parameter,  P)  may  be  obtained  from  an  accurate  zero  energy  trial  function  it 
can  be  said  that  the  benefits  obtained  from  doing  the  zero  energy  calculation 
based  on  a  minimum  principle  may  be  utilized  in  the  case  of  slightly  higher 
energies  as  well.  Numerical  calculations  for  the  problem  of  electron  scattering 
by  atomic  hydrogen  (for  P  and  D,  as  well  as  S-vave  scattering)  are  now  being 
performed.  I*  is  hoped  that  with  the  aid  of  the  present  techniques  the 
discrepencies  regarding  the  total  low  energy  e"*H  scattering  cross  section,  which 
have  recently  been  reported,  may  be  resolved. 
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Table  I.  Results  of  variational  calculations  of  the  singlet  e~H  scattering 
length,  in  units  of  the  Bohr  radius.  In  the  first  column  the  values  obtained 
by  Borowitz  and  Greenberg,  who  used  what  we  have  called  the  'inappropriate' 
normalization  [see  Eq.  (2.I4)],  are  given.  The  values  which  appear  in  the 
second  column  were  obtained  by  the  present  authors  by  converting  the  Borowitz- 
Greenberg  data  to  the  Kohn-Hulthen  farm  [  'appropriate'  normalization;  see 
Eqs.  (2.2)],  thereby  obtaining  a  rigorous  upper  bound.  The  fact  that  a  slight 
improvement  is  gained  in  this  conversion  is  in  agreement  with  the  discussion 
in  the  text.  It  was  possible  to  obtain  a  bound  because  the  u.  used  by  the 
above  authors  included  a  bound  state  component  multiplied  by  a  variational 
parameter. 

NORMALIZATION 

'INAPPROPRIATE'      'APPROPRIATE' 

NO  POLARIZATION  8. 16  8.1b 

APPROXIMATION 

POLARIZATION  7.75  7.63 

CONSIIERED 
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